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چکیده 

در حال حاضر باقلا به‌عنوان یکی از محصولات مهم در منطقه خاورمیانه. افريقاء چین و نقاطی از اروپا و استرالیا در 
تغذیه انسان و دام مورد توجه قرار دارد. مدل. ابزاری است که ما را در تفسیر و درک دنیایی که در آن زندگی می‌کنیم. 
پاری می‌کند. پیش‌بینی دقیق فنولوژی گیاهان زراعی از ویژگی‌های ضروری مدل‌های شبیه‌سازی به‌شمار می‌رود. از 
مدل‌ها می‌توان در بهبود مدیریت تولید گیاهان زراعی برای پیش‌بینی تاریخ‌های احتمالی برداشت یا پیش‌بینی عملکرد 
نهایی. يا به‌صورت فعال‌تر» برای پیش‌بینی مطمئن زمان وقوع حوادث فنولوژیکی استفاده نمود. مدل استفاده‌شده در این 
مطالعه. مدل ارائه‌شده توسط سلطانی و همکاران برای نخود بود. برای پارامترهای ورودی مدل. از مقادیر به‌دست‌آمده از 
ازمایش‌های مزرعه‌ای در دو سال زراعی ۱۳۹۰-۹۱ و ۱۳۹۱-۹۲ که در مزرعه تحقیقاتی دانشکده علوم کشاورزی و منابع 
طبیعی دانشگاه گنبد کاووس در شرایط مطلوب دیم و عاری از افات» علف‌های هرز و بیماری‌های باقلا اجرا گردیده بود. 
استفاده شد. ۲۴سناریوی مختلف شامل ترکیبات مختلفی از تراکم و تاریخ کاشت در دو سال زراعی برای مدل تعریف 
گردید. نتایج نشان داد که مدل قادر است مراحل فنولوژیک گیاه باقلا را به‌جُز روز تا سبزشدن به‌خوبی پیش‌بینی کند. این 
پیش‌بینی در مراحل منتهی به برداشت باقلا از بالاترین دقت برخوردار بود. به‌نحوی که بالاترین ضریب تبیین با مقدار ۰/٩۱‏ 
مربوط به تعداد روز تا برداشت باقلا بود. مطالعه نشان داد که مدل فوق نمی‌تواند به‌خوبی حداکتر سطح برگ. عملکرد. 
عملکرد بیولوژیک و شاخص برداشت را در شرایط آب و هوایی گنبد و تحت شرایط آزمایش‌های مزرعه‌ای پیش‌بینی نماید. 


واژه‌های کلیدی: توصیف مدل. حبوبات شبیه‌سازی» مدل 


مقدمه 

باقلا با نام علمی م 1004 ۷/25/4 متعلق به خانواده 
6 می‌باشد. با احتمال زیاد منشأً باقلا جنوب‌غربی آسیا 
گزارش شده است (2005 ,۳1۵۵56181 ۷]2[0008). باقلا با 
دارابودن ۲/۴ ۲درصد پروتئین در دانه خشک باعثت گردیده تا 
درحال‌حاضر به‌عنوان یکی از محصولات مهم در منطقه 
تقاهرسانهه افر تابن وتقاطی از ارفیاه انلیا قز تقد یه آشمان 
و دام مورد توجه قرار گیرد (101010 2002 ,. 61). سطح 
زیرکشت باقلا در جهان و ایران به‌ترتیب معادل ٩/۲میلیون‏ 
و۳۰هزار هکتار می‌باشد. مناطق عمده کشت باقلا در ایران. 
استان‌های گلستان. خوزستان مازندران و گیلان می‌باشد 
(2006 ,۳۲0956181 ۷]2[2000). سطح کشت باقلا در استان 
گلستان در سال زراعی ٩۳‏ - ۰۱۳۹۳ ۲۲۸۸هکتاره میزان تولید 
آن ۲۰۷۶۱تن( ۲۵درصد تولید باقلا در کشور) و عملکرد غلاف 


*نویسنده مسئول: ۵۱۷۵600.60۳0 تصعطهتتله 


۳۸ 


سبز آن ۰۷۴٩کیلوگرم‏ در هکتار گزارش شده است. باقلا گیاه 
مناسبی در تناوب با غلات است. زیرا موجب شکستن چرخه 
بسیاری از بیماری‌های متداول غلات می‌شود و جمعیت نماتدها 
را کاهش می‌دهد (2017 ,.۸1 1 212130060. 

مدل, ابزاری است که ما را در تفسیر و درک دنیایی که در 
آن زندگی می‌کنيم. یاری می‌ کند. دانشمندان و مهندسان از ائواع 
مدل‌ها به‌عنوان ابزاری برای درک پدیده‌های مورد مطالعه استفاده 
رک مد اش ممادنه با وهای از اد لانق است ک ۸ 
رفتار هر سیستم را به‌طور کمّی توصیف می کند ( ,.61 2 ۳0016 
#۶ صصتایم ۷2۵ :1996 مصمصصعناهو عک تتقاممله :1996 
3 .). پسیش‌بینی دقیق فنولوژی گیاهان زراعی از 
ویژگی‌های ضروری مدل‌های شبیه‌سازی به‌شمار می‌رود 
(2004 ,10۵01 :2006 .]۵ ۶۶ ]601 

تولید و تسهیم ماده خشک در مدل‌های شبیه‌سازی گیاهان 
زراعی تا حد زیادی به وسیله زمان‌بندی مراحل نمو تنظیم 
می‌شود ( 6 50[]201 : 2005 متصفااهک ک 12۵10هک1 تحعطف۴! 
6 ./0) درک چگونگی واکنش نموّ فنولوژیک باقلا به طول 
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روز و دماء ایجاد مدلی برای پیش‌بینی زمان وقوع مراحل مختلف 
نموّ این گیاه زراعی را میسر می‌سازد. انجام آزمایش‌ها و تحقیقات 
علمی بر روی محصول زراعی و یافتن شرایط مطلوب رشد آن‌ها 
مستلزم هزینه‌های گزاف و صرف وقت و دقت زیادی می‌باشد و 
این موضوع برای کشورهای درحال‌توسعه به علت کمبود امکانات 
و پژوهشگر دشوار است. امروزه با استفاده از مدل‌های زراعی و 
امکانات گسترده‌ای که در زمینه نرم‌افزارهای کامپیوتری به‌وجود 
آمده است. می‌توان با پیداکردن روابط میان فرایندهای رشد و 
یافتن شرایط مطلوب رشد به نتایج قابل‌توجهی دست یافت و 
آن‌ها را برای شرایط دیگر تعمیم داد و وضع جدید را پیش‌بینی 
نمود. مدل‌ها محدودیت‌های جغرافیایی و محیطی را کاهش داده 
و در مورد گیاهان زراعی برای ارقام مختلف قابل تعمیم می‌باشند 
(1996 ,صمصجع‌تاهگ ک عتقاممتو :1996 باه ۶۶ 6ا00ظ. 

استفاده از مدل‌ها بسیار سریع بوده و با درنظرگرفتن گیاه 
زراعی به عنوان یک سیستم پویاء امکان پیش‌بینی تغییرات این 
سیستم نسبت به زمان را فراهم می‌آورند. علاوه بر آن, از این 
مدل‌ها می‌توان برای معرفی گیاهان جدید در یک منطقه نیز 
استفاده کرد. پیش‌بینی برای وضعیت تولید مواد غذایی در یک 
منطقه یا جهان و اتخاذ سیاست‌های لازم با توجه به تغییرات 
محیطی نیز از دیگر مزایای مدل‌سازی است ( ,۵ 2 ۳0066 
6 ,56112120 6 تتقاممتگ :1996). 

از مدل‌ها می‌توان در بهبود مدیریت تولید گیاهان زراعی 
برای پیش‌بینی تاریخ‌های احتمالی برداشت یا پیش‌بینی عملکرد 
نهایی, یا به‌صورت فعال‌ت برای پیش‌بینی مطمتن زمان وقوع 
مراحل فنولوژی به نحوی که کود. تنظیم‌کننده‌های رشد. علف 
کش‌ها و کنترل بیماری‌ها در مناسب‌ترین زمان به‌کارگرفته شوند. 
استفاده نمود. بسیاری از مدل‌سازان از مدل‌های خود به‌عنوان 
وسیله‌ای برای ارزیابی خطرات موجود در تولید استفاده کرده‌اند 
(1999 ,۱۷۲۷۵۵00 6 1101217). این امر با بررسی واکنش 
عملکرد (مدل) به آمار درازمدت اقلیمی یک منطقه صورت 
می‌گیرد و در نتيجه آن بهترین تاریخ کاشت. تراکم» فاصله ردیف؛ 
رقم زراعی و غیره تعیین می‌شود. از مدل‌ها می‌توان در بررسی 
کمّی اثر ویژگی‌های زراعی بر روی رشد و عملکرد گیاهان در 
محیط‌های خاص استفاده کرد. در مورد مسائل به‌ن زادی گیاهان 
زراعی همچون دورگ‌گیری دیررسی و زودرسی بهبود دانه‌بستن؛ 
مورفولوژی گل‌هاء افزایش رشد تک‌دانهها و زمان گلدهی می‌توان 
2 مفل‌ها استفاده کرق کعیین تاتسن عملکرف ععطقه‌ای: کسک 
به مدیریت آبیاری و ارزیابی اثرات تغییر اقلیم از کاربردهای دیگر 
مدل‌ها می‌باشند ( ,107۵0 :1996 مصفحصعناع5 ک عتق01ط1 
9 ۱۷۲۱۵۱۵۱0۲ 6 51801211 :2004 
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(2017) ,۵1 ۰۶ له ررظ در مطالعه‌ای که روی لوبیاچیتی 
انجام دادند. نتیجه گرفتند که می‌توانند از مدل ترابی و سلطانی 
برای پیش‌بینی عملکرد و مراحل فنولوژیکی این گیاه استفاده 
کنند و بیان داشتند که این مدل می‌تواند به خوبی این پارامترها 
را بیان کند. 

استفاده از یک مدل شبیه‌سازی برای پیش‌بینی رشد و 
عملکرد. به پارامترهای ورودی آن مدل در منطقه موردنظر 
بستگی دارد. مدل‌های شبیه‌سازی گیاهان مثل سایر ابزارها دارای 
نقاط قوت و محدودیت‌های خاص خود هستند. بنابراین» حصول 
نتیجه مطلوب از به‌کارگیری مدل‌ها مستلزم شناخت مناسب از 
توانایی‌ها و محدودیت‌های آن‌ها می‌باشد. هدف از این تحقیق. 
ایجاد یک مدل ساده برای پیش‌بینی مراحل فنولوژیک, حداکثر 
تتظیم یگ سگرن و شاخمی رعاش اقا بت تفت ابط 
آب‌وهوایی گنبد با استفاده از مدل نخود بود. 


مواد و روش‌ها 
آزمایش‌های مزرعه‌ای 

این آزمایش در دو بسال فراعی ۱۳۹۰-۹۱ و ۱۳۹۱-۹۲ کار 
مزرعه تحقیقاتی دانشکده علوم کشاورزی و منابع طبیعی دانشگاه 
گنبد کاووس در شرایط مطلوب دیم و عاری از آفات» علف‌های‌هرز 
و بیماری‌های باقلا اجرا گردید. گنبدکاووس در ۵۵درجه و 
۲دقیقه طول جغرافی‌ایی و ۲۷درجه و ۱۶دقیقه عرض 
جغرافیایی» در ارتفاع ۴۵متر بالاتر از سطح دربا قرار دارد. میزان 
بارش سالانه حدود ۰ ۰ ۵میلی‌متر است. قبل از اجرای آزمایش از 
عمق صفر تا ۲۰سانتی‌متری خاک نمونه‌برداری انجام شد و 
خصوصیات فیزیکوشیمیایی خاک تعیین شد (جدول ۱) که بر 
اسان تنج حاصاه باق اک لوم برس سیافی نود 

آزمایش به‌صورت فاکتوریل در قالب طرح بلوک‌های کامل 
تصادفی با چهار تکرار و در دو سال انجام شد. تراکم کاشت در 
چهار سطح (۸ ۰۱۲ ۱۶ و ۲۰ بوته در مترمربع) و تاریخ کاشت در 
سه سطح (۱۲آبان ۰۱۳۹۱ ۱۲آذر ۰۱۳۹۱ ۱۳دی ۱۳۹۱ و 
۷آذر ۰۱۳۹۰ ۲۴آذر ۱۳۹۰ و ۱۷دی۱۳۹۰) بود. برای ایجاه 
تراکم‌های مختلف. فاصله بین ردیف‌ها ۰ ۵سانتی‌متر و فاصله روی 
ردیف‌ها ۵ ۲سانتی‌متر برای تراکم ۸بوته در مترمرسع. 
۷سا مر براق فراک ۲ ابوتق در متزمریی: ۱۲۵ سای سر 
برای تراکم ۱۶بوته در مترمربع و ۱سانتی‌متر برای تراکم ۰ آبوته 
در مترمربع درنظرگرفته شد. میزان کود توصیه‌شده در زمان 
کاشت. ۰ ۲کیلوگرم سوپرفسفات تریپل و ۵۰کیلوگرم کود اوره در 
هکتار بود. بذور قبل از کاشت با سم کربوکسی تیرام به میزان دو 
در هزارنه ضدعفونی شده و در ردیف‌های ایجادشده در عمق 
۵سانتی‌متری کشت شدند. چون این آزمايش در شرایط مطلوب 


راحمی کاریزکی و حسینی؛ مدل‌سازی مراحل فنولوژیک.../ پژوهش‌های حبوبات ایران/ سال(دوره)۱۱ شمارة!. نیمة اول ۱۳۹۹ صفحة ۳۸-۴۸ 


دیم و عاری از آفات» علف‌های‌هرز و بیماری‌های باقلا انجام شد. 
هر دو تا سه هفته یک‌بار برای مبارزه با آفات ازجمله شته سیاه 
باقلا (۳۵0۵6 ۸۵0۳15) از سم ایمیداکلوپرامید (کنفیدور) به 
میزان ۲۵۰سی‌سی در هکتار و به‌منظور جلوگیری از بیماری‌های 
قارچی مانند لکه‌شکلاتی (۳۵0۵6 15ا0۱ظ) و برق‌زدگی 
(م۳۵0 ۵2۹۵00]۵) از سم بنومیل به نسبت دو در هزار 
استفاده گردید. به‌خاطر انجام عملیات تهیه زمین مناسب. قبل از 


از آن به‌صورت هفتگی عملیات وجین به‌صورت دستی انجام شد. 


پس از سبزشدن بذور در مرحله ۴ تا ۶برگی تنک‌کردن انجام شد. 
ضقانت اندارهگیر‌شهه هر این کی شاب سای مراخلن مه 
فنولوژیک ازجمله روز تا گلدهی (161» روز تا غلاف‌دهی (153), 
روز تا دانه‌بندی (155» روز تا شروع رسیدگی (57) و روز تا 
رسیدگی کامل (158) بود که در فصل رشد روی ؛ ابوته مشخص 
بر اساس روش (1977) 022۷16085 6 76107 انجام شد. به این 
نحو که ۱۰بوته انتخاب و با یک نخ قرمز مشخص گردید و سپس 
با فاصله هر دو روز یک‌بار مراحل فنولوژی مورد بررسی قرار 
گرفت. 


جدول ۱- مشخصات خاک محل اجرای آزمایش (عمق ۳۰۰-۰ سانتی‌متری) 
(سن 30 ۲۵0 0 ط)۵ع0) ععهام اماصمعصتهوی صا ممتامننهاه فطل لز50 :120161 


مقدار مشخصه 
ص۷۵ عناعنها هط 

173 گوگرد (قسمت در میلیون) (00۳) 1000۲ 

11 هدایت الکتریکی (دسی زیمینس) (0) 

16 اسیدیته ([۲۳) 

105 درصد مواد خنثی‌شونده ۷۵ ..] 

116 کربن آلی (درصد) (۶0) 00ات عنصحع:0۵ 

03 نیتروژن (قسمت در میلیون) (0۳۲0) 0عع0تانل( 

123 فسفر قابل‌جذب (قسمت در میلیون) (ظع) ۲مطحردمطم داعاتام۸ 
414 پتاسیم قابل‌جذب (قسمت در میلیون) (005) اندعداهم زاعانامدطاض۸ 
30 رس (درصد) (36) نود 

62 لای (درصد) (66) 116 

8 ماسه (درصد) (6۵) 200 


توصیف مدل 

مدل استفاده‌شده در این مطالعه مشابه با مدل نخود 
2۵ عراعم [1-]99[۷ توسط (1999) ,۵ ۶ تطدااه5 
بود. در این مدل جهت پیش‌بینی پارامترهای مختلف فنولوژیکی 
نیاز به مجموعه‌ای از اطلاعات پایه است که در جدول ۲ آورده شده 
نت برای استفاده از این مدل. ابتدا اقدام به تعریف سناریوهایی بر 
اساس تاریخ کاشت و تراکم مورد مطالعه در آزمایش مزرعه‌ای 
گردید و در نهایت خروجی مدل (خروجی‌های مربوط به مراحل 
شاخص برداشت) با داده‌های واقعی به‌دست‌آمده از آزمایش 
مزرعه‌ای مورد ارزیابی قرار گرفت. 

داده‌های هواشناسی و اقلیمی بر اساس داده‌های مربوط به سال 
های ۱۳۲۹۰ تا ۱۲۹۲ گنبد مورد استفاده قرار گرفت. برای اجرای مدل. 
۴ ۲سناریو شامل ترکیبات مختلف تاریخ کاشت و تراکم مورد استفاده 
در آزمایش‌های مزرعه‌ای تعریف گردید. جهت مقایسه داده‌های شبیه 
سازی‌شده با داده‌های به‌دست‌آمده از آزمایش‌های مزرعه‌ای» از 
شاخص‌های ارزیابی. ضریب تبیین (87 جذر میانگین مربعات خطا 


( 12197 جذر میانگین مربعات خطای نرمال‌شده ( 01۳۹۳ و 
خط ۱:۱ استفاده گردید. هرجه مقدار **1 حاصل از رگرسیون خطی 
توابع و خط ۱:۱ به یک نزدیک‌تر باشد» نشان‌دهنده همبستگی بالا بین 
مقادیر شبیه‌سازی‌شده و مشاهده‌شده است و نشانه توصیف بهتر مدل 
در شبیه‌سازی صفت موردنظر می‌باشد. علاوه بر آن مدل ایده‌آل آن 
است که در رگرسیون مقادیر شبیه‌سازی‌شده در مقابل مقادیر 


اندازه گیری‌شده ضرایب ۰ و 02۱ به دست آیند. 


نتایج حاصل از مدل بیانگر پاسخگویی خوب مدل برای 
مراحل مختلف فنولوژیکی گیاه باقلا در شرایط گنبد بود. مدل به 
غیر از مرحله روز از کاشت تا سبزشدن, سایر مراحل «روز از کاشست 
تا گل‌دهی. روز از کاشت تا غلاف‌دهیء روز از کاشت تا دانه‌دهی. 
روز از کاشت تا رسیدگی فیزیولوژیک و روز از کاشت تا برداشت) را 
به‌خوبی برآورد کرد (جدول ۲). 


1. ۵0۲ ۳820 50۲۵۲۵ 0۲ 
2. ۱0۳۳۲۵117260 1۲001 12620 ۹0۲۵۲۵ ۲ 
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جدول ۲- پارامترهای ورودی موردنیاز مدل 
[۱۱۵06 0۲ عمامه 2 اامرصز ۱۱۵۵060 2۰ 12016 
واحد پارامتر (مخفف) 
اتصرنا ۱ 


داده‌های هواشناسی و اقلیمی 
ععل عافصناه 0صه ۱۷۸/626۲ 


درجه سلسیوس(0؟) 5تاژوام) 168766 حداکثر دمای روزانه (۲2)076)1[۷]۸26ه جع ولندل صیاصمتنه]۱۷ 
درجه سلسیوس (0؟) تالواع) 168766 حداقل دمای روزانه (1(۷]1۳) 6تنااهتهمججها راتمل حصتصتصن۷( 
مگاژول در مترمربع در روز 02۷ ۳6۲ 24 تعاعحه و50۲2 وم عانامز2ع۱۷]6 تابش خورشیدی روزانه «بلهکاع) عمن20ن۲20 تا50 
میلی‌متر 1۷۲1111۳06167 بارندگی روزانه ((۵ الهکصنه۲ «رانددز 
ویژگی‌های خاک 
دمته 220 9011 
ندارد 0[ آلبیدو (ظ/6۸1) 40عحلض۸ 
میلی متر بر میلی ۳ ضریب زهکشی (0۳3۸۵1(۲) ۲ماععا معفمنه:ظ۳ 
۲مامصتاا ۱۷۲ مه تمامصصتااز۷ 
سانتی‌متر بر سانتی‌متر کسر حجمی‌رطوبت در نقطه پژمردگی دنا 
تمعن تهج مامصتامون اتصنا همینا 0عصنحل اه اصماجم مه امد عتتاعصتاآن ۷ 
سانتی‌متر بر سانتی‌متر کسر حجمی‌رطوبت در نقطه ظرفیت زراعی «ألل) 
امن تهج مامصصتامون اتصنا م1۱ من اه اصماصم تمه تمد عتتاعصصتان ۲۷ 
سانتی‌متر بر سانتی‌متر کسر حجمی‌رطوبت در نقطه اشباع (15۸1) 
امن تهج مامصتامون ۵ 2 اجماومه مه تمد عتتاعصصتاه ۲۷ 
ندارد 0[ شماره منحنی خاک (()) اناد وت 
مدیریت زراعی 
اصعصصم‌عمصهحه متصمط0عخ 
روز سال از اول ژانویه 12۳11277 ]15 عطا مه کهع۷ ۵۶ ردرز[ تاریخ کاشت (۳۲۵۷) 216 507102 
بوته در مترمربع ۳616 ۹00۵۲6 06۲ )۳۱ تراکم ((ظ) 0510۲ع([ 
میلی‌متر رطوبت خاک در موقع کاشت (15۷۷7) 
تمامصتا۷ 9 آ ۱۱۱۱۶۱۰۱ 
ندارد کسر آب قابل‌تعرق خاک برای انجام آبیاری 101/۷1) 
۱۵ [۷۵ 0۱4ععط) ممناهع 1۳1 


جدول ۳- پارامترهای ریاضی و آماری مربوط به صفات مختلف مورد مطالعه به‌دستآمده از مدل و آزمایش‌های مزرعه‌ای 


کاجه‌صتهمون 610 0صه م۵۵ ما وا 0مز0ناای مات عبامره۲۱ عظ) م6 0ماها۲۵ فهمامصصونوم اف‌تاعتاهاه 20 )۱۷۲2/۵۵/1 :3 12016 


پارامترهای آماری میانگین حداقل حداکثر 
مهو آه‌تاعتاهاد ۱/۹۹ نصا تا ۱/۵۹۵۱ صفت 
۵ ۵ 1۳۱ ۳ پیش بینی‌شده واقعی پیش بینی شده واقعی پیش بینی‌ شده واقعی ۳2۳۵ 
او ۱ 002 ۴ ۸02 ۴( ۸21 
029 623 019 9 36 22 5 14 43 دقرب ری 
6 10 121۹ 
ز تا گلد 
012 وو2 089 108 96 88 87 131 111 0 
110002 0) ورزر[ 
ی 
018 ۹47 031 121 114 100 101 141 131 ۱ 
۷ 660 0۲ تاو ما ورد[ 
۳۷ 
016 51 034 127 118 108 103 147 134 بر روز ۶ دنه 
۵ 5660 0 12215 
تا 9 ۵ ۰ ۳ 
012 584 09 158 150 138 128 115 118 روز تا رسید کی فیزیولوژیکی 
۲امص آهععمامتووطام م) ویرچرز 
ز تا برداشت 
011 623 091 166 157 136 133 183 190 4 
۷۲ 1۵ 1215 
029 8 ۰ 04 ۱58 253 082 ۳ 27 3271 سطح برگ ماکزیمم 
(جص) مهد ک2ع1 حصتحصت‌جه]۱۷ 
0۳ ۱ 
0 .۰ 3242 0017 ...1902 3253 1294 2025 298 5194 ی وت 
(حصلع) ۷۱614 660 
شاخ اشت 
029 53 0006 501 463 44 34 56 522 4 
(6۵) 1006۶ ۳۲۱۵۲۷65۲ 
یک حّ 
014 همم 009 381 4903 252 305 497 728 عملکرد بیولوژیک (گرم در مترمریج) 


حصع) 10عزر اعتعمامظ 


۳۱ 
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روز از کاشت تا گل‌دهی: جهت بررسی این پارامتر روز 
تاگل‌دهی مشاهده‌شده از مزرعه باروز تاگل‌دصی 
پیش بینی شده توسط مدل مورد ارزیابی قرار گرفت. مقادیر 
مشاهده‌شده در آزمایش مزرعه‌ای با حداکثر و حداقل 
معادل ۱۱۱ و ۸۷روز و ميانگین ۹۶روز بود. در صورتی که این 
مقادیر در نتایج پیش‌بینی‌شده معادل ۱۳۱ و ۸۸روز و میانگین 
۱۰۸ بود. مقایسه بین روز تا گلدهی واقعی 9 پیش‌بینی‌شده 
نشان داد که داده‌ها در دامنه ۱۵درصد از خط ۱:۱ قرار دارند 
بین این مقادیر در خصوص روز تا گلدهی برابر ۰/۸۳۹ به‌دست 
امك: جذر میانگین مربعات خطا 9 مقدار نرمال‌شده ق به تر تیب 
معادل ۳/۹۹ 9 ۲ روز بود (جدول ۲). 

روز از کاشت تا غلاف‌دهی: نتایج نشان داد که مدل. روز 
تا غلاف‌دهی را نیز به خوبی پیش‌بینی کرد. به‌نحوی که ضریب 
خطای نرمال‌شده. بین داده‌های به‌دست‌آمده از آزمایش‌های 
مزرعه‌ای و داده‌های پسیش‌بینی‌شده از مدل به‌ترتیسب 
معادل ۰۰/۸۱ ۲/۹۹ و ۰/۱۲ بود (جدول۲). مقایسه داده‌های 
حاصل مقادیر اندازه گیری‌شده و پیش ‌بینی‌شده نسبت به خط 
۱ در محدوده ۱۵درصد این خط قرار گرفتند (شکل ۱) که 
آن هم بیانگر پیش‌بینی خوب مدل از خصوصیت فنولوژیک روز 
تا غلاف‌دهی باقلا بود. 
حاصل از آنالیز رگرسیون خطی توابع بین مقادیر روز تا شروع 
رشد دانه اندازه گیری‌شده و شبیه‌سازی‌شده گیاه باقلا در 
شرایط گنبد نشان داد که مدل. ارزیابی خوبی از ایین پارامتر 
داشت (۰/۸۴- **). علاوه بر آن» جذر ميانگین مربعات خطا و 
دما کین رات حف ای تفه بوربه تیب با 
مقادیری معادل ۵/۱ و ۰/۱۶ بیانگر پیش‌بینی خوب این پارامتر 


اندازه گیری‌شده روز تا شروع رشد دانه و مقادیر پیش‌بینی‌شده 
ی نیز نشان داد که تمامی‌داده‌ها در محدوده ۱۵درصد خط 
۱ قرار داشتند (شکل ۱). 


روز از کاشت تا رسیدگی فیزیولوژیکی: مقادیر روز تا 
رسیدگی فیزیولوژیک به‌دست‌آمده از مزرعه در یک محدوده‌ای 
از مقادیر ۱۲۸ تا ۱۷۸روز و میانگین ۱۵۰روز به‌دست آمده در 
مقابل نیز اعذاه به‌دست‌آمده از مدل نیز تقریبا فر دامته‌ای 
مشابه با اعداد واقعی به‌دست آمد. به‌نحوی که حداکثر و 
حداقل‌روز تا رسیدگی فیزیولوژیک به‌دست‌آم ده از مدل 


۳۲ 


بت یسب ای ۱۷۵ ۱۸ ارو اتکی ارو تیک 
آمد. در خصوص این پارامتر ضریب تبیین» جذر میانگین 
مسا ها و خر مالکيق ترسات خط ای تفه 
به ترتیب برابر ۰۰/۹ ۵/۸۴ و ۰/۱۲ به‌دست آمد که همگی نشان 
دهنده پیش‌بینی خوب مدل می‌باشد (جدول۳). مقادیر روز تا 
رسیدگی فیزیولوژیک واقعی و پیش‌بینی‌شده مدل نیز در 
محدوده ۱۵درصد از خط ۱:۱ قرار دارند (شکل ۱). 

روز تا برداشت: نتایج مطالعه نشان داد که روز تا برداشت 
پارامترهای مورد مطالعه برخوردار بود. به‌نحوی‌که مقدار 
حداکثر روز تا برداشت که توسط مدل پیش‌بینی‌شده بود. 
معادل ۱۸۳روز و حداقل این عدد معادل ۱۳۶روز و با میانگین 
۶ روز پیش‌بینی گردید. در مقابل حداکثر» حداقل و میانگین 
روز تا برداشت به‌دست‌آمده از آزمایش‌های مزرعه‌ای به‌ترتیب 
نوی بر فک مسبت اسف 
حاصل آنالیز رگرسیون خطی توابع بین روز تا برداشت 
اندازه گیری‌شده و مقادیر پیش‌بینی‌شده برابر ۰/۹۱ بود که این 
مقدار به خوبی نشان می‌دهد که مدل. مقادیر مناسبی را از 
این پارامتر پیش‌بینی نموده است (جدول ۲). علاوه بر آن. سایر 
تااش هی امن کمنا تط ۵ که ابید تام مطالبت 
فوق می‌باشد (شکل ۱). 

(1997) 1620۲ ک 100021) به منظور پیش‌بینی 
عملکرد گندم در دشت‌های آبی منطقه پنجاب هندوستان با 
استفاده از مدل و داده‌های اقلیمی پنج‌ساله. عملکرد دانه. ماده 
خشک کل, مراحل فنولوژیک. گلدهی و رسیدگی فیزیولوژیک 
گندم را شبیه‌سازی کردند. مقایسه نتایج شبیه‌سازی‌شده با 
نتایج مشاهده‌شده در مزرعه نشان داد که مدل به طور رضایت 
بسن عهاکت ۵و هم کش کم اقفر یتک 


به دست‌آمده حاصل تابع رگرسیون بین مقادیر داده‌های 
پیش‌بینی‌شده و واقعی نیز به‌دلیل کمتربودن از ۰/۶ بیانگر این 
مطلب بود (۰/۴- 02 (جدول ۲. (2008) ,۵1 6 266 نیز 
بیان کردند که مدل مورداستفاده در مطالعه ایشان قادر نبود 
تا به خوبی حداکثر شاخص سطح برگ را پیش‌بینی کند. 
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شکل ۱- مراحل مختلف فنولوژیکی مشاهده‌شده در مقابل مراحل مختلف فنولوژیکی پیش‌بینی‌شده 
خل جویزه هرگ سل ۸ و خطاوط تنظه‌شین بالا و جانیی یه خر قیب بانط ۱۵ مات ش باق 


قمع اهمزعهامصمصه م0 ۳6010۲60 کنام۲۵ ممعه)ر اهماممامصمصاص م06 0۵۵۹۵۰۲۷۵۵ :1 ۲12۰ 
۰ ر1:)0.85 200 1:1.15 هتم دما عطقم عع10۷ هه عممرجرنا صرح 1:1 فصن و فصن فیامتصتاومت) 
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عملکرد دانه: مقایسه آموزه‌های آماری مربوط به مقایسه 
بین عملکرد به‌دست‌آمده از آزمایش‌های مزرعه‌ای و عملکرد 
پیش بینی‌شده در مدل نشان داد که مدل تحت شرایط مطالعه 
قادر نیست تا پیش‌بینی خوبی را از عملکرد به ما بدهد. 
به‌عبارتی ضریب تبیین» جذر میانگین مربعات خطاو جذر 
میانگن مربعات خطای نرمال‌شده به‌ترتیب برابر ۰۰/۰۱۷ ۲۱/۴۲ 
و ۰/۱ به‌دست آمد (جدول۳). علاوه بر این. در خصوص نمودار 


۰.0 ۰ 500 ۰ ۰.400 ۰ ۰ 300 ۰ ۰ 200 ۰ 100 0 
خی ۲۱6۱0 ادسا۸ 


0 700 600 500 400 300 200 100 0 
نع ۲6۱0 ای‌نعمامنط ادسا‌ظ۸ 


(مصای ۵۱ ۲۱۳۵۵۵0 


(خصصی) ۷۱۵۱0 اد‌زعمامنها ۲۱۳۲۵۵۲۵۲۵۵۵ 


بین عملکرد واقعی و عملکرد پیش‌بینی‌شده نیز مشخص گردید 
که داده‌ها تقریبا خارج از محدوده‌ی ۱۵درصد خط ۱:۱ قرار 
دارند (شکل ۲). نتایج به‌دست‌آمده از اين مطالعه با یافته‌های 
(2015) 107201 > 3221671 در خصوص باقلا برای استان 
گلستان مطابقت نداشت. (2017) ,47 6 ۸۲۷۵112 نیز در مورد 
لوبیاچیتی به نتایج مشابه دست پیدا کردند. 
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شکل ۲- حداکثر سطح برگ (۰)۸ عملکرد (). شاخص برداشت (6) وعملکرد بیولوژیک () مشاهده شده در مقایل مقادیر 


پیش‌بینی‌شده آن‌ها 
خط توپر, معرف خط ۱:۱ و خطوط نقطه‌چین بالا و پایین به ترتیب بیانگر خط ۱:۱/۱۵ و ۰/۸۵: ۱ می‌باشد. 
کاصتامحظ فتاه صمط) میعه۲ ۲۱۵۱0 امزممامزه 20 0) 1026 امم2۲۲۷ ور۵) ۲۱6۱0 ررض) همه آهم۱ تفه م۳۲۵0 2۰ ۲۱2۰ 


۰ ,1:0.85 هه 1:1.15 هه معصناً 0عطفقل م1027 هه تممروربا هه 1:1 فصن ور فصن فیامتصتاجمت) 
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راحمی کاریزکی و حسینی؛ مدل‌سازی مراحل فنولوژیک.../ پژوهش‌های حبوبات ایران/ سال(دوره)۱۱ شمارة!. نیمة اول ۱۳۹۹ صفحة ۳۸-۴۸ 


شاخص برداشت: نتایج نشان داد که مدلء پیش‌بینی 
مناسبی را از شاخص برداشت باقلا برای شرایط گنبد نمی‌دهد. 
با وجود قرارگرفتن بیشتر داده‌ها در محدوده ۱۵درصد خط 
۱ پارامترهای آماری بیانگر عدم توجیه مناسب مقادیر 
پیش ‌بینی‌شده شاخص برداشت توسط مدل بود. به‌طور مثال. 
ضریب تبیین برابر ۰/۰۰۶ به‌دست آمد (جدول ۲ و شکل ۲). 

عملکرد بیولوژیک: عملکرد بیولوژیک به‌دست‌آمده در 
آزمایش مزرعه‌ای در دامنه‌ای از حداقل معادل ۲۰۵ گرم در 
مترمربع تا حداکثر ۷۸۷گرم در مترمربع و میانگین ۴۹۰/۳ گرم 
در مترمربع به‌دست آمد. مدل نیز مقادیر حداقل» حداکثر و 
مبانگین عماک ره بیولوزیک را به‌ترقیب مساال 3۵۳ ۳۹۷ و 
۱ گرم در مترمربع محاسبه نمود. مقایسه بین این دو سری 
داده نیز نشان داد که مدل در خصوص این پارامتر نیز همانند 
عملکرد اقتصادی قادر نیست تا پیش‌بینی مناسبی را ارائه 
نماید. پارامترهای آماری شامل ضریب تبیین» جذر میانگین 
مربعات خطا و جذر میانگین مربعات خطای نرمال‌شده نیز 
بیانگر این مطلب بود (جدول ۲). 


نتیجه گیری 

نتایج حاصل از پیش‌بینی مدل نشان داد که مدل می‌تواند 
به خوبی مراحل مختلف فنولوژیکی گیاه را برای شرایط آب و 
خی گنبم بیس پیت هید پیت بیتی سول در مرح 
فنولوژیکی منتهی به زمان برداشت به‌ خصوص پیش‌بینی‌روز تا 
رسیدگی فیزیولوژیک وروز تا برداشت از دقت خیلی بالاتری 
نسبت به سار مراحل فنولوژیکی برخوردار بود (جدول ۱). در 
مقابل نتایج این مطالعه حاکی از عدم پیش‌پینی مناسب 
پارامترهای عملکردی شامل عملکرد اقتصادی. شاخص برداشت 
و عملکرد بیولوژیک گیاه باقلا در شرایط مطالعه بود. علاوه بر 
آن» مدل پیش‌بینی مناسبی را از حداکثر سطح برگ به ما ارائه 
نکرد. شاید بتوان دلیل اين امر را اين‌گونه بیان نمود که چون ما 
از مدل نخود جهت مدل‌سازی برای گیاه باقلا استفاده نمودیم و 
همچنین بعضی از پارامترهای موردنیاز در مدل به‌دلیل عدم 
وجود این پارامترها برای گیاه بقلاءناگزیر از پارامترهای نخود 
استفاده گردید. این امر باعث شد تا مدل نتواند به خوبی بعضی 
از پارامترها را توجیه کند. 
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6 1:1 عطا م0 20 ۵ 1570 ۵۶ فصو عط مر مج وملج۷* عوعط ۶ه له ع2ط ٩۳00۱760‏ 061مصظ مه مصم 0عصتهاط0اه 
00610160 مملامصهامنه مط) اقطا م9 .عصتااعو ۵00 60 ره 0۳6016060 ۷۷۵11 2150 2۵061 معط )62 5۳00۱760 و16تاوع۲ م1" 
0 ۲0۳ ۵0021260 اه عطا ممع۵ مه هنود صحعص )۲۵0 ۵۵1260 4ص تمه عتقتا۹0 .39828 ]۲۵0 
۵۵8 ۲6۵5060۷۵۱۷۰ ,12.0 20 ,2.99 ,0.81 60 لهیاهع ع7۷۵ [2۵06 عط ۵۶ 0202 0۲6010060 عط مه فاصممهمرم 
۸ 217 0۳۲۵0۲ 50۷۵۲۵ 10820 )۲۵0۵ بطمتاقصهامونه ۵۶ اصمل‌گلهمه مطا رفظ آهمتعمامتو روط فطا م6 رقل عطا 
6 0۶ ه0۲6010010 20۵00 ۵ اصعوع:6۵: ۵1 من ,۲۵۹۵6۵۵0۷۵ ,12:0 2۳00 ,5.84 ,0.9 ۷۷۵۲۵ 6170۲ 50۲۵۶۵ 102621 )۲۵۵ 


حصمی.ممطه ۵ تصمطه‌نله تتمطاسم مصنل‌صموم‌ممن* 
۳۷ 


۳۱۵۱6۵ 0۴۲ افصسول صهنصه] 
(۸۲۵۱۵ 6862۲0۱ع1) 38-48 .0 و2020 و1 .۱0 و11 ,۷۵01 


مصمتاعص عط که ممتومع‌وع: عفعمناً عمط عصترواقمه وه 0عصنقادان موه رح تصعئ‌تگلمی ممتتقمهاوه 16 ,۱9۸۵061 
0۰ مه اهنا0ع ومتالد۷ 010)60عن عطا مج )2۷65 ما ره ممع۲7عه 


جمتعصل‌صم) 

مفطا از بصهه‌ها 1202 مط ۵۶ قمعه)ه اهم1ع۱0مصمطاج فط) )۵۲6016 ۷۷۵11 ۱0نامي 390061 فطل عقط) 07۵0و فاناوع1 فط 1 
صجهه م6 ۱62018 ومع2)ه مطا صا [۱6۷۵ ]وعطعنظ مط) 21 فه۷۵ ممنا6010ع0 ونط1 .,عمممهمصصه عطا مه رد مط) ۵۶ ممتاممعرم 
5 0۶۴ ناه فط 10۶ 0صنام مه 0.91 ۵۶ متام و طز ممتاهمهامونه ۵۶ )مملهمتالهمن )وعطعنط مط کقطا 50 راوع۷ت2ظ 
,1306 2762 ک2ع۱ مرتاصرل همه مط) ۴۲۵ )016ع۲ج مصصوی [ع۵0 عماج فط فص ۵0 مه تاه نان باوعتقط موه هه 
۰ 27۱ نت0 رع نا مه فصمت)تصمن تمطهع ۳۷ مصبال مد »رع1۳0 و2۷۵ قمع ۷1۵10 101021621 ,۷1610 


مم نصا مممتامتع وم ۱۷۱۵0۵1 بام۱۷]۵0 رفصفه‌ظ :1۵۲۱۷۵۲۵۵ 


۴۳۸ 


